Multiple shape coexistence in the nucleus ^^Zr 
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Abstract 

We study the low- lying energy spectrum of the rp-process waiting point nucleus ^^Zr with state-of-the-art beyond mean 
field methods with the Gogny DIS interaction. Symmetry restoration and configuration mixing of axial and triaxial 
shapes are included in the calculations. Five 0+ states corresponding to different nuclear shapes are obtained below 2.25 
Q^MeV and several rotational and 7- bands built upon them are identified. Nevertheless, these states do not modify the 
^ /3-decay half-life having a negligible effect in the rp-process. A good agreement with the available experimental data is 
obtained. 
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The study of shape coexistence in nuclear systems is 
one of the most active topics of current research. In par- 
ticular, shape coexistence have been observed and/or the- 
oretically predicted in the proton-rich region with mass 
number A ^ SO [1-8]. These nuclei are characterized by 
having at least one 0+ excited state in the low lying en- 
ergy spectrum that can be associated to a shape different 
from the ground state one. In this region, such an effect is 
mainly due to the competition between the closure of the 
pf subshells and the intrusion of some levels of the gg /2 and 
6^5/2 shells whenever the deformation is increased. Both 
schematic and self-consistent interactions show this partic- 
ular behavior [9-12]. As an example, we display in Fig. 1 
the single particle energies calculated with the Gogny DIS 
interaction [13] as a function of the quadrupole deforma- 
tion ^2- The Fermi level crosses the N = Z = 40 spherical 
gap and, increasing the deformation, the pi/2 and oi"- 
bits are depopulated while ^9/2 and ^5/2 levels enter below 
the Fermi energy. This structure favors the appearance of 
different minima in the potential energy surface and the 
competition and/or coexistence of different shapes at low 
energies. The nucleus ^^Zr {N = 40 protons and Z = 40 
neutrons) belongs to this region of potential shape coexis- 
tence. 

Experimentally, there are only six levels measured in 
this isotope, being identified as members of the ground 
state rotational band [14, 15]. From the theoretical point 
of view, this isotope has been previously studied using 
microscopic- macroscopic models [9-11], Monte Carlo shell 
model [16] and self-consistent mean field approximations 
with Skyrme [12, 17, 18], Relativistic [19, 20] and Gogny 



interactions [21]. There are also calculations with these 
interactions including beyond mean field effects, like axial 
angular momentum projection [22] with additional mixing 
of axial shapes [23], parity projection [24] and in the frame- 
work of the collective Hamiltonian [25]. Most of these 
works are focused on the ground state properties of ^^Zr 
and not on its excitation spectrum. For the description 
of shape coexistence several ingredients are required such 
as the triaxial shapes with the corresponding angular mo- 
mentum projection and shape mixing [26], and, addition- 
ally, the particle number projection [27]. The latter is im- 
portant: 1) to remove the possible spurious contributions 
coming from eigenstates of the particle number operators 
with different eigenvalues from the desired one, and 2) if 
the projection is performed before the variation [28, 29] 
to avoid spurious phase transitions in the weak pairing 
regime ( pairing collapse ) and to provide the optimal en- 
ergy since one applies the variational principle to a wave 
function with sharp particle number at variance with the 
HFB case. 

Recently, we have presented [30] a beyond mean field 
approach, more precisely a Symmetry Conserving Configu- 
ration Mixing (SCCM) theory, which considers linear com- 
bination of quadrupole deformed wave functions projected 
to particle and triaxial angular momentum quantum num- 
bers. The finite range density dependent Gogny interac- 
tion was used and a model application was done to study 
the nucleus ^^Mg with a small configuration space. We 
now apply this state-of-the-art SCCM theory for the first 
time to a case requiring a large configuration space, namely 
the study of the prospective axial and triaxial shape co- 
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Figure 1: (color online) Single particle energies for neutrons (protons 
follow a similar pattern with energies about 11 MeV higher) as a 
function of the quadrupole deformation ^2 calculated with the Gogny 
DIS interaction. The Fermi level is represented by a thick dotted 
line. 



existence in the nucleus ^^Zr. The high predictability of 
the Gogny interaction provides an added bonus to the cal- 
culations. As it will be shown below, the approach used 
here improves the one given in Ref. [25] although the un- 
derlying interaction is the same. There, particle number 
projection is not performed and the overlaps needed for 
the shape mixing are assumed to be gaussian to build the 
collective hamiltonian. These approaches are valid under 
certain conditions but not in the whole triaxial plane. 

Proton-rich nuclei are also of paramount relevance to 
describe the synthesis of elements in the so-called rapid- 
proton capture process (rp-process) which occurs in X-ray 
bursts produced in binary systems composed by a massive 
star and a neutron star [31, 32]. Along the rp-process path 
there are some even-even N = Z nuclei such as ^^Ge, ^^Se 
and ^^Kr with a small proton capture cross-section and 
a relatively long beta decay half-life which are known as 
waiting points. Hence, final abundances of the elements 
produced in these bursts and also the time-scale of the 
nucleosynthesis process depend crucially on the structure 
of these specific nuclei. The nucleus ^^Zr is one of these 
rp-process waiting points although the beta decay half-life 
of its ground state is shorter than in the cases referred 
above [33]. However, the appearance of isomeric states in 
the low-lying energy spectrum could open additional /3+ 
decay channels modifying the effective half-life of this nu- 
cleus, speeding-up or delaying the rp-process. 



The nuclear shape is a semiclassical concept which can 
be associated with a mean field description. Therefore 
the optimal approach to describe nuclear shapes is a self- 
consistent mean field theory (Hartree-Fock-Bogoliubov, HFB, 
in our case). By constraining the pertinent operators (the 
quadrupole mass operators Q2fi) one is able to produce 
different shapes at will. In our case the well known defor- 
mation parameters (/32,7) are defined as: 
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being tq = 1.2 fm and A the mass number. 

HFB wave functions (w.f.) -I^^st)" eigenstates 
of symmetry operators like the angular momentum (AM) 
and the particle number (PN) and a restoration of such 
symmetries is mandatory. For finite systems the PN sym- 
metry must be imposed from the very beginning and the 
w.f.'s 1^/327) are found by solving the equations resulting 
of a variation after the PN projection (PN-VAP approx- 
imation) [34], i.e., minimizing the projected energy with 
the mentioned constraints. With these w.f.'s one is able to 
calculate PN and AM projected potential energy surfaces 
(PNAMP approach): 
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where /, M, K are the AM and the z component of the vec- 
tor / in the laboratory and intrinsic frame respectively and 
pi,N,z ^ pi^pNpZ and Pjo (Pmk in general) and P^(^) 
are the AM and the neutron (proton) projectors [35]. 
The energy eigenstates are not in general eigenstates of 
the shape operators and a mixing is expected if the local 
minima are close in energy. Therefore our final approach 
is a SCCM calculation (with shape and K mixing) in the 
framework of the Generator Coordinate Method (GCM): 
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<j = 1, 2, ... labels the states for a given value of the angular 
momentum ((7 = 1 corresponds to the Yrast states). The 
coefficients /x/327'^ ^ are found by solving the Hill- 

Wheeler- Griffin equations that are deduced from the mini- 
mization of the energy calculated with the non-orthogonal 
set of wave functions {Puk^^ P^\^ ^2^)} [35]. The so- 
lution of these equations provides the energy levels and 
all the information needed for calculating, for example, re- 
duced transition probabilities, the so-called collective wave 
functions (square of the probability density, normalized to 
1, of finding the state (/, a) with given deformation param- 
eters (/32,7)), decomposition of the states in etc. More 
details about how the observables are computed within 
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Figure 2: (color online) Potential Energy Surfaces -normalized to their absolute minima respectively- calculated for ^°Zr using the Gogny 
DIS interaction: (a) only particle number symmetry is restored -before variation- and (b) the rotational symmetry is also restored (/ = 0), 
i.e., the PNAMP approach. Spatial densities for each minimum of (b) are also shown. Contour lines are separated by 0.2 MeV (dashed) and 
1 MeV (continuous). The color scale is the same for (a) and (b). 



this framework can be found in Reference [30] (and ref- 
erences therein). In this work, HFB states \^f32-f) con- 
serve parity, time-reversal and simplex symmetries and 
they are expanded in a spherical harmonic oscillator ba- 
sis which includes nine major oscillator shells. Concern- 
ing the input parameters chosen in these calculations, we 
use A^GCM = 60 intrinsic states to form a triangular mesh 
in the (/32,7) plane. Finally, although the exploration of 
the octupole degree of freedom is beyond the scope of the 
present letter, we have checked the influence of {3^ per- 
forming axial HFB calculations in the range of quadrupole 
deformations studied here. As a result, we do not find any 
solution where the parity symmetry is broken. 

To study the possible shape coexistence in this nu- 
cleus we start with analyzing the potential energy sur- 
face (PES) in the triaxial (/32,7) plane calculated with 
the PN-VAP approximation. In Fig. 2(a) we observe a 
very rich structure with five minima located at deforma- 
tions (/32,7) = (0,0°), (0.53,0°), (0.37,20°), (0.47,32°) 
and (0.2,60°). The axial minima (7 = 0°,60°) can be di- 
rectly related to the gaps observed in the single-particle 
energies represented in Fig 1. In this approach, the abso- 
lute minimum corresponds to the spherical point which is 
in contradiction to the rotational behavior observed exper- 
imentally. The other minima are, respectively, at 4.0, 4.4, 
5.7 and 4.1 MeV above the spherical one. Furthermore, 
this energy landscape displays a rather gamma-soft be- 
havior. The barriers between the prolate, the triaxial and 
the lowest oblate minimum are smaller than 1 MeV. Nev- 
ertheless, this structure is expected to change whenever 
beyond mean field effects, including angular momentum 
projection and shape mixing, are taken into account. 
In Fig. 2(b) the PES in the PNAMP approach (/ = 0) 
is shown. We also represent the spatial density distribu- 
tion of the intrinsic states that correspond to the different 



minima found in the surface to give a better idea about 
their shapes. We observe that, except for the spherical 
point, the rotational correction induced by the symmetry 
restoration produces a significant lowering of the energy 
in the whole plane. The prolate and oblate minima in the 
PN-VAP approach have been displaced to triaxial shapes 
and slightly larger j32 deformations. In fact, the absolute 
minimum corresponds now to a moderately triaxial de- 
formed shape (/32,7) = (0.55,7°) while the spherical one 
has been displaced to a slightly prolate shape, being now 
at ~ 1.4 MeV above the absolute minimum. This result 
goes in the right direction to reproduce the experimen- 
tal data. The effect is somewhat similar to the ^yVg2Q 
case in the island of inversion where the deformed con- 
figuration becomes the ground state only after restoring 
the rotational symmetry [36]. Nevertheless, we have to in- 
clude the possible infiuence of shape and K mixing in the 
calculation because the barriers between neighboring tri- 
axial minima remain of the order of 1 MeV or less. Hence, 
we perform SCCM calculations including particle number 
and angular momentum projection (Eq. 3) to obtain the 
energy levels and electric quadrupole reduced transition 
probabilities B{E2). The predicted spectrum is shown in 
Fig. 3. The states are sorted in bands according to the val- 
ues of the highest transition probabilities B{E2) between 
those levels. Highest values of B{E2) within the different 
bands are also displayed in Fig. 3 while several interband 
transitions are given in Table 1. Additionally, the distri- 
bution of probability of having a given deformation (/32, 7) 
for the band-head states is also plotted on the right in 
Fig. 3. This distribution corresponds to the square of the 
collective wave functions in the triaxial plane. Here we 
observe very well defined maxima showing that the con- 
cept of shape is still valid even after the shape mixing. 
The energy ordering of these maxima coincides with the 
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Figure 3: (color online) Left: Experimental [14, 15] and theoretical spectra for ^°Zr. B{E2) values are given in e'^b'^ Right: Square of the 
collective wave functions for the band heads. 



minima of Fig. 2(b) with the exception of the spherical 
and the oblate ones which have been exchanged. Interest- 
ingly, the shape fluctuations bring back the "almost" pro- 
late and oblate minima to pure prolate and oblate shapes 
as in Fig. 2(a). In addition, the distribution of probability 
remains practically the same for the rest of the members 
of a particular band and corroborates the stability of the 
different minima. The distribution in K for the states be- 
longing to a band is also similar as it is shown in Table 2. 
Furthermore, the most intense electric quadrupole transi- 
tions, both between intraband and interband states, are 
related to similar distribution in K of the initial and final 
states (see Tables 1, 2 and Fig. 3). 

In Fig. 3, on the left, we observe four 0+ excited states 
below 2.25 MeV, built on top of different shapes, coex- 
isting with the ground state. The ground state band is 
predicted to be a rotational band built on top of an ax- 
ial prolate deformed configuration with the usual sequence 
/ = 0,2,4,6,... (g.s. in Figs. 3 and 2(b)). The collective 
wave functions have their maxima at (/32,7) = (0.55,0°) 
that corresponds to the minimum found in Fig. 2. Com- 
paring these results with the experimental energies we ob- 
tain a rather good agreement. This shows again the rele- 
vance of restoring the rotational invariance in this nucleus 
where a spherical ground state is found at the PN-VAP 
approach. The theoretical values of the energies tend to 
follow a more rotational behavior /(/ + 1) than the data 
and the value for the B{E2^ 2i ^ Oi) seems to be overes- 
timated according to the trend of the experimental data 
along the isotonic = 40 chain [37]. These effects could 
be related to the fact that we are not including time- 



reversal symmetry breaking intrinsic states. The angu- 
lar momentum projection of cranked states would com- 
press the spectrum -affecting only the excitation energies 
for states with 1^0- and would modify also the value 
of the transition probabilities [38]. All the states belong- 
ing to the g.s. band correspond to rather pure K = 
configurations (see Table 2). According to a semiclassical 
picture [35], a 7- vibrational band with K = 2 is expected 
to be present with an / = 2+ state as the band head and 
AI = 1. These states are found indeed at 2.4 MeV as we 
can see in the figure (7-g.s. in Fig. 3, see also Table 2). 
The transition probabilities along this band are large and, 
since i^-mixing is very weak, those to the g.s. band very 
small. 

The first excited band built on top of the O2" state {Eq+ = 
0.76 MeV) has also a rotational character based on a tri- 
axial configuration (1st exc. in Figs. 3 and 2(b)). The 
distribution of probability is concentrated at (/32,7) ~ 
(0.4,20°). As we see in Table 2, this is mostly a K = 
band although some mixing with K = 2 components is 
present, especially in the / = 62", 83" states. We can also 
find a quasi-7 band partner with the 2^ state at 1.32 
MeV as the band head, and with the probability peaked 
at the same deformation as its partner (7- 1st exc. and 
Fig. 3). The even / states of this band have a considerable 
= 0, 2 mixing while the odd / levels are rather pure 
K = 2 states. Correspondingly, we obtain large B{E2) 
values for A/ = 2 transitions between the levels of this 
band and rather strong inter-band transitions to the first 
excited band. 

We also find another band and its quasi-7 band partner 
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Table 1: Selected theoretical interband electric quadrupole reduced 
transition probabilities B(E2). Following Fig. 3 the short notation 
gsb stands for ground state band, leb for 1st excited band, lej for 
1st excited 7-band, etc. 



la 


K = 


\K\=2 


\K\=4 




K = 


\K\=2 


|K| =4 


Oi 


1.000 






26 


0.022 


0.978 




2i 


0.997 


0.003 




33 


0.000 


1.000 




4i 


0.995 


0.004 


0.001 


48 


0.005 


0.980 


0.015 


61 


0.990 


0.010 


0.000 


54 


0.000 


0.976 


0.024 


O2 


1.000 






23 


0.210 


0.790 




22 


0.889 


0.111 




3i 


0.000 


1.000 




42 


0.804 


0.191 


0.005 


43 


0.352 


0.613 


0.035 


62 


0.725 


0.260 


0.014 


5i 


0.000 


0.946 


0.054 


82 


0.638 


0.326 


0.034 


64 


0.201 


0.767 


0.024 


03 


1.000 






25 


0.267 


0.733 




24 


0.707 


0.293 




32 


0.000 


1.000 




44 


0.462 


0.498 


0.040 


45 


0.239 


0.312 


0.449 


63 


0.570 


0.369 


0.058 


52 


0.000 


0.770 


0.230 


83 


0.525 


0.383 


0.084 


65 


0.262 


0.069 


0.581 


05 


1.000 














27 


0.777 


0.223 












49 


0.808 


0.095 


0.097 











Table 2: Decomposition of the collective states appearing in Figure 3 
in their < 4 components 



with strong K mixing built on top of the triaxial minimum 
of the PES at (/32,7) = (0.5,30°) (2nd exc. and 7-2nd exc. 
in Figs. 3 and 2(b)). The excitation energies of the band 
heads are E{0^) = 1.4 and £^(2+) = 1.8 MeV respectively. 
In the quasi-7 band the largest values for the B{E2) cor- 
respond to A/ = 2 and also A/ = 1 transitions with an 
odd / initial state. The states of these two bands are also 
strongly connected to the first excited quasi-7 band. Fi- 
nally, we obtain two more 0^ states located on top of the 
spherical and oblate minima of the PES at excitation en- 
ergies of 1.5 and 2.24 MeV, respectively. 

From the astrophysical point of view, the appearance 
of isomeric states in the low-lying energy spectrum could 
modify the effective decay half-life towards the ^^Y 
isotope. However, we obtain half-lives of the different ex- 
cited 0+ states much shorter than the typical half-lives for 
a /3+ decay process. Furthermore, the excitation energies 
of these states are too high to allow for their population 
at the typical temperatures of the rp-process T ~ 1 GK 
(equivalent to E - 100 keV) [32]. Therefore, the /3+ de- 
cay of ^^Zr is expected to occur from its ground state only 
which has an experimental half-life of ~ 4.1 s [33]. 

In summary, we have applied state-of-the-art beyond 
mean field methods to study the structure of the rp-process 
waiting point nucleus ^^Zr. We find five 0^ states below 
2.25 MeV that correspond to different minima in the triax- 
ial potential energy landscape which constitutes a unique 
example of multiple shape coexistence in the nuclear chart. 
This feature only occurs if full triaxial angular momen- 



tum projection is performed. However, we found that due 
to their short lifetimes the effect of the excited 0+ states 
on the beta-decay of ^^Zr and therefore on the rp-process 
is negligible. An experimental observation of the excited 
0+ states predicted by our calculations, once radioactive 
beam intensities will be high enough to perform for exam- 
ple Coulomb excitation studies of ^^Zr, will certainly be of 
highest interest. 
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